Interrelation of work function and surface stability: the case of BaAl 4 
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The relationship between the work function ($) and the surface stability of compounds is, to 
our knowledge, unknown, but very important for applications such as organic light-emitting diodes. 
This relation is studied using first-principles calculations on various surfaces of BaALt. The most 
stable surface [Ba terminated (001)] has the lowest <E> (1.95 eV), which is lower than that of any 
elemental metal including Ba. Adding barium to this surface neither increases its stability nor 
lowers its work function. BaALt is also strongly bound. These results run counter to the common 
perception that stability and a low <£> are incompatible. Furthermore, a large anisotropy and a stable 
low-work-function surface are predicted for intermetallic compounds with polar surfaces. 
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I. Introduction 

The long-standing problem of the precise relationship 
between the work function (<I>) and the stability of met- 
als has become pressing with the increased application 
of electron-emitting materials in technology. Two exam- 
ples are vacuum electronic devices like cathode-ray tubes 
(CRTs) and cathodes for organic light-emitting diodes 
(OLEDs). In CRTs, a thin layer of a low-<!> metal is 
often present on top of a cathode made from a struc- 
turally stable material to enhance its electron-emitting 
properties^ Electron injection into OLEDs strongly de- 
pends on the cathode work function,- while the lifetime 
of the device can be limited by the stability of the cath- 
ode material. OLEDs with an alloy interface between the 
cathode and the polymer have been found to be superior, 
in terms of lifetime and luminosity, to those with single- 
element cathodesA^ Thus the relationship between $ 
and stability^ is crucial. However, it is poorly under- 
stood, especially for more complex metals. 




FIG. 1: Schematic illustration of Smoluchowski smoothing of 
the electron charge density (striped area) at an open surface. 
The grey area is formed by bulk Wigner-Seitz cells. Signs 
indicate the net charge. The dipole thus formed lowers the 
work function. 

The general rule for elements is that a low work func- 
tion and high stability are incompatible: the element 
with the lowest work function, cesium (<I> = 2.14 eV)f& is 
highly reactive and has a low melting temperature. No- 
ble metals (silver/gold/platinum) on the other hand are 
hardest to oxidize, but their $ is at least twice as large 
(4.25/5.1/5.65 eV)£ It is generally believed that this 
must always be the case, i.e., that a low $ implies loosely 



bound electrons that easily mediate reactions. The work 
function surface anisotropy, however, can be quite large: 
for tungsten, it is of the order of 1 eVi The anisotropy 
was already theoretically addressed by Smoluchowski. 8 
According to his model, at a "more open"— surface, re- 
laxation of the electrons "smoothes" the surface charge. 
A dipole moment is built up that lowers <&>iS This is 
schematically shown in Fig. ^ Although the work func- 
tion is lowered for a more-open surface, the surface energy 
increases and stable, low-work-function surfaces are for- 
bidden by the model. Hardly any extensions to the model 
have been suggested sincejii certainly not for more com- 
plex metals. Experimental results for transition metal 
alloys suggest the so-called alloy effect^ which implies 
that the stability and <& of alloys interpolate between 
those of the constituting elements. As a consequence, 
it is believed that stable, low-work-function surfaces are 
not possible for more complex metals either. 

Nevertheless, the work function and stability of 
compounds mainly constitute a terra incognita. In 
this paper, the relationship between <I> and stability of 
compounds is studied for the first time by calculating 
ab initio the structural relaxation, work function and 
surface energy of various surfaces of BaALi. The mo- 
tivation to study the BaAl4 system was the successful 
use of Ba-Al-alloy cathodes in OLEDs and the high 
melting point of this compound (> 1000° C). In the 
Ba-Al phase diagram^ the melting point of BaALi is 
the highest, much higher than those of the constituents, 
which indicates a strongly bound structure. Moreover, 
we show that the work function for one of the crystal 
surfaces is very low. 

II. Ab Initio Calculations 

The BaALi crystal structure^ is depicted in Fig. [3 It 
is body-centered tetragonal, with alternately three alu- 
minum layers (1A1, 2A1 and 3A1) and one barium layer 
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(IBa) in the [001] direction. Both Al and Ba atoms are lo- 
cated approximately according to their elemental crystal 
structures (bcc and fee, respectively). The (100) surface 
is stoichiometric and there are four (001) surfaces. They 
are constructed by cutting the bulk above the accordingly 
labeled layers in Fig. [3 Two other (001) surfaces are also 
considered: "2Ba" (one bcc barium layer added to IBa) 
and "iBa" (half a barium layer removed from IBa). 




FIG. 2: Body cen- 
tered tetragonal unit 
cell of BaAU contain- 
ing 2 formula units. 
The labels refer to 
the layer and atom 
type. 1A1 and 3A1 
atoms are equivalent 
in the bulk. 



From bulk calculations we find that a = 4.56 A, that 
c = 11.39 A , and that the remaining free parameter in 
the structure, the height of the third aluminum layer, 
is 0.381 • c. These values compare very well (deviations 
w 1%) with those obtained experimentally^ The bulk 
density of states (DOS) shows a quasi-gap just below the 
Fermi level, in accordance with the previously calculated 
DOS of BaAU and other isoelectronic compounds with 
the same structured The valence electrons are mainly 
located on Al atoms. The binding energy (£hi n d) i s 
1.42 eV/F.U. with respect to the elemental bulk metals. 
Together with the quasi-gap, it hints at BaALi's stability. 

The first-principles calculations were carried out using 
density functional theory (DFT) in the local density 
approximation (LDA)A£ii£ with generalized gradient 
corrections (GGA)*i£ We used the total energy and 
molecular dynamics program called VASP (Vienna 
Ab-initio Simulation Package) jASiiSl which has the 
projector-augmented-wave method (PAW)22i2I imple- 
mented. Nonlinear core corrections^ were applied for 
both barium and aluminum. A semi-core of Ba 5s and 
bp electrons was included. The Kohn-Sham orbitals 
were expanded in plane waves with cutoffs of 18 Ry. 
The Brillouin zones for the (100) and (001) surface 
calculations were sampled with 1x16x8 and 12x12x1 
Monkhorst-Pack 2 ^ k point grids, respectively, resulting 
in eight and 21 k points, respectively, in their irreducible 
parts. The (periodically repeated) unit cell for the (100) 
surface calculation contained a slab with a thickness of 
seven bulk unit cells and 16 A vacuum. The supcrcells 
for the (001) surface calculations contained slabs with 
thicknesses of six formula units (F.U.) and 20 A vacuum. 



III. Work Function 
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FIG. 3: The BaAl 4 (001) IBa surface. Laterally averaged 
charge density (p solid line, arbitrary units) and (electro- 
static) potential (V dashed line, eV relative to Ef) as function 
of position perpendicular to the surface (z). The positions of 
Ba and Al layers are indicated with light and dark spheres, 
respectively. At —5.12 eV is the average bulk potential and 
at 2.79 eV is the highest potential in the bulki— The work 
function is 1.95 eV. 

The work function is defined as the amount of energy it 
takes to extract electrons from a metal, i.e., bring them 
from the Fermi level to the vacuum. At locations that 
are microscopically far from the material, but macroscop- 
ically near it, <I> is surface dependent. The work func- 
tion at large distance is then an average over the various 
surfaces. 25 The work function of a surface is calculated by 
constructing a superccll with a slab of material, with only 
this surface, and vacuum. Fig.[3]iflustrates this using the 
(001) Ba terminated surface. About 10 A empty space 
suffices for the electrostatic potential 2 ^ to converge to its 
vacuum value (Vv ac ). A thick slab, however, is required 
for the Fermi level (Ep) to be accurate. A more efficient 
method is to link the surface calculation to a bulk one 
where Ep is very accurate. The average potential in the 
bulk (< V >buik) is set equal to the average potential in 
the middle of the slab. Accuracies of a few hundredths 
of an eV can thus be achieved with only six to eight bulk 
unit cells, depending on direction^ The work function 
of the IBa surface is 1.95 eV. 

The work functions of the other BaAU surfaces are cal- 
culated analogously. That of the (100) surface is reduced 
by structural relaxation from 3.54 eV to 3.43 eV. The 
work functions of the (001) 1A1, 2A1 and 3A1 surfaces 
are 4.38 eV, 4.48 eV and 3.98 eV, respectively. Relax- 
ation has negligible effect. The work functions of the 
various surfaces are summarized in Fig. 0] As plotted 
there, covering the Al surface with half a monolayer of 
barium (thus forming the ^Ba surface) already greatly 
reduces the work function (2.29 eV). The IBa work func- 
tion is reduced even more to a surprisingly low value of 
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FIG. 4: Work functions (eV) of the BaAl 4 (001) surfaces. 
Experimental bulk values for barium and aluminum are indi- 
cated at the borders. Lines connect the data points. $ for 
the (100) surface before and after relaxation is inserted. (The 
dotted line is VJnax)" 



1.95 eV. This is below the $ of any element^ An ad- 
ditional layer of barium (thus forming the 2Ba surface) 
increases $ (2.13 eV) again. Experimental bulk values for 
Al (4.2 eV) and Ba (2.32 eV) are plotted for referenced 
These compare favorably with calculated work functions 
of bcc-barium [2.36 eV (100) and 2.27 eV (111)] and fcc- 
aluminum [4.34 eV (100) and 4.17 eV (111)]. 

The huge (> 2 eV) variation in the (001) surface work 
function can be (qualitatively) explained from the dif- 
ference in the atomic electronegativities. Ba is less elec- 
tronegative (0.9)2£ than Al (1.5) j2£ effecting charge trans- 
fer from barium to aluminum. The resulting surface 
dipole moment decreases the work function of the Ba 
surface to even under the elemental barium bulk value 
and increases $ of the Al surfaces considerably. The de- 
pendence of $ on Al coverage of the IBa surface can be 
understood if one realizes that the 1A1 and 3A1 surfaces 
are more open than the 2A1 surface. According to the 
Smoluchowski model, <& must then be lowered. The de- 
crease in $ during relaxation of the (100) surface can also 
be understood by combining charge transfer with Smolu- 
chowski smoothing. The barium atoms at the surface 
move out (« 0.1 A), as they favor an environment with 
a lower charge density, while the aluminum atoms at the 
surface tend to get closer together, as they favor an en- 
vironment of a high charge density. We conclude that a 
mono-layer coverage yields an extreme work function. 
IV. Surface Energy 

The (relative) stability of a surface is determined 
by the difference in surface energy (7) between two 
surfaces^ The surface energy is calculated as the energy 
of a slab with only this surface minus the energy of the 
bulk equivalent (formed by merging the [periodically re- 
peated] slabs together), divided by two times the surface 
area, because a slab has two sides. This method was used 



for the (100) and (001) ^Ba surfaces. Slabs for the other 
(001) surfaces are non-stoichiometric and so there is no 
equivalent bulk. Calculating 7 now requires reservoirs of 
Ba and Al. Assuming thermodynamic equilibrium, the 
chemical potential of the aluminum (fiAi) and that of the 
barium (/iBa) reservoir are linked to the total energy per 
F.U. of BaAk (£ bulk ): 

-Etmlk = 4 • HAl + MBa- (1) 

The energy of, for example, the IBa surface, as a function 
of //Ba, follows from the total energy of a slab with these 
surfaces exclusively (-EiBa)i^ 

7 1Ba (W3a) = [ElBa ~ #F.U. • £ b ulk - MBa]/2A S , (2) 

where A$ is the surface area. 
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FIG. 5: Surface energies (eV/nm ) for the BaAU surfaces as 
function of the barium chemical potential (eV). The poten- 
tial ranges from the value for Al bulk to that of Ba bulk. 
The difference is the binding energy of BaAU (equation has 
been used). For the (100) surface the effect of relaxation is 
indicated. The (001) IBa surface is the most stable. 

The surface energies are drawn in Fig. Structural 
relaxation of the (100) surface lowers its energy from 
4.96 eV/nm 2 to 4.82 eV/nm 2 . The energy of the (001) 
iBa surface equals 13.3 eV/nm 2 and is outside the range 
of the plot. For the other (001) surface, 7 depends on 
the barium chemical potential. In the figure, /ZBa varies 
over the thermodynamically allowed range between the 
Ba and Al bulk chemical potentials— as phase separa- 
tion occurs outside this range. The single barium surface 
is the most stable one in the entire region and so other 
(001) surfaces cannot be formed. 

The unusual stability of the low-<!> surface is explained 
in the same way as the relaxation of the (100) surface. 
Since barium is the less electronegative element, it 
favors a low electron density and makes the most stable 
surface. Additional barium layers at the surface decrease 
this stability, as Ba-Ba bonds are less strong than Al-Ba 
bonds, especially when the barium in contact with 
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aluminum is (partially) ionized. 

V. Conclusions. 

As both the low $ of the IBa surface and its stability 
followed from Ba's lower electronegativity we come to 
a remarkable prediction: For an intcrmetallic compound 
with polar surfaces, like BaAl4, the most stable surface 
has the lowest work function and relaxation can only de- 
crease it further. 

To summarize, we used first-principles calculations on 
various surfaces of BaAl4 to study the interrelation of 
work function and stability. The bulk work function of 
BaAl 4 is 2.79 eV^i its anisotropy (1.5 eV) is huge, and 
the most stable surface (Ba terminated 001) has the low- 
est $ (1.95 eV), which is even lower than that of pure 
Ba. Adding barium to this surface neither increases its 
stability nor lowers its work function. The binding en- 
ergy of 1.42 eV/F.U., the quasi-gap in the DOS, and its 
melting temperature of over 1000 °C indicate BaALj's 



stability. These results contradict the common percep- 
tion that stability and a low $ are incompatible. They 
also run counter to the effect that alloy work functions 
are in between those of the constituting elements. Its sta- 
ble, very-low-$ surface and the stable structure probably 
make BaAU a good electron-emitting material. Further- 
more, a stable low-work-function surface promises to be 
general for intermetallic compounds with polar surfaces. 
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